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SECTION  I 


INTRODUCTION 

The  interaction  of  a  transverse  jet  with  an  oncoming  stream 
is  a  fundamental  and  important  aerodynamic  phenomenon  in  many 
flow  regimes.  In  the  subsonic  domain,  examples  include  the 
discharge  of  gases  from  smokestacks,  mixing  and  combustion,  and 
VTOL/STOL  technology.  For  supersonic  flows,  the  characteris¬ 
tics  of  the  transverse  jet  interaction  are  important  to  control 
and  maneuverability  of  aerospace  vehicles  as  well  as  supersonic 
combustion.  However,  our  knowledge  of  this  complex  physical 
phenomenon  is  still  very  limited  and  urgently  needs  continued 
research  efforts. 

Investigators  have  considered  many  aspects  of  the  inter¬ 
acting  transverse  jet  problem,  including  jet  dynamics 
(References  1-9),  entrainment  and  mixing  (References  1-3,  6, 
10),  induced  pressure  distributions  (References  8,  11,  12), 
shock  interaction  (References  8,  11,  13-15),  and  numerical 
solutions  (References  5-7,  13,  16).  Few  investigators  have 

attempted  to  calculate  the  complete  three-dimensional  flow 
field.  Even  fewer  have  attempted  numerical  simulations 
utilizing  the  ensembled  compressible  Navier-Stokes  equations. 
None ,  to  our  know  ledge ,  have  studied  the  nature  of  the  flow 
field  when  the  jet-issuing  body  is  rotating  with  respect  to  its 
principal  axis . 

For  the  case  of  a  supersonic  primary  stream  and  a 
non-rotating  body,  the  flow  field  is  characterized  by  the 


inviscid-viscous  interaction  similar  to  that  induced  by  a 
protuberance  (References  17,  18).  Complex  patterns  of  shock 
interaction  together  with  three-dimensional  vortical  structures 
are  created  as  the  flow  is  diverted  around  the  jet.  The 
phenomenon  of  interference  which  forms  separation  regions  is 
well  documented  for  many  types  of  flow  field  obstructions 
(References  17,  18) .  In  the  case  of  a  jet  in  a  hypersonic 
primary  stream,  the  domain  of  influence  is  significantly 
reduced.  Therefore,  the  gradients  of  flow  properties  are  much 
steeper  and  the  scales  of  the  resultant  flow  structures  are 
changed.  The  three-dimensional  separated  flow  must,  however, 
obey  topological  rules  which  maintain  the  continuity  of  the 
flow  field  (References  17,  19-21).  In  order  to  capture  the 

detailed  flow  field  structure,  a  large  amount  of  computer 
resources  is  required.  When  the  jet  issuing  body  rotates  about 
its  longitudinal  axis,  there  is  no  plane  of  symmetry  to  reduce 
the  computational  domain,  and  the  grid  resolution  problem 
becomes  more  severe.  The  first  portion  of  the  present  study 
establishes  the  criterion  for  numerical  rest'  ition  of  the 
three-dimensional  interaction  which  is  generated  by  a  super¬ 
sonic  jet  issuing  from  a  blunt-nosed  ogive  cylinder  into  a 
hypersonic  stream. 

The  major  portion  of  the  present  study  concerns  the  cal¬ 
culation  and  correlation  of  the  flow  field  past  the  jet  issuing 
cylinder  rotating  about  its  longitudinal  axis  under  identical 
test  conditions.  The  numerical  simulation  is  not  dynamically 
coupled;  the  ogive  cylinder  is  ass  ’med  to  remain  at  zero  angle 


of  attack  at  all  times.  Two  approaches  are  possible  in 
describing  the  rotational  fluid  motion.  The  first  method  is  to 
define  the  coordinates  in  the  Newtonian  frame  and  generate  the 
instantaneous  grid  system  according  to  the  motion.  The  other 
approach  is  to  describe  the  fluid  motion  in  the  rotating  frame 
of  Reference  22.  This  latter  approach  is  more  suitable  to  our 
present  purpose.  First,  calculation  in  the  rotating  frame  of 
reference  is  simpler  because  the  grid  need  not  be  rotated  at 
each  time  step.  Second,  since  the  deformation  tensor  is 
invariant  under  rotational  transformation,  calculation  in  the 
rotating  frame  of  reference  still  allows  us  to  determine  vis¬ 
cous  effects  (such  as  the  shear  flow  pattern)  while  maintain¬ 
ing  an  easily  implemented  computer  code.  Finally,  the  addi¬ 
tional  acceleration  due  to  the  non-rectilinear  motion  can  be 
singled  out  easily  to  assess  the  impact  on  the  complete  flow 
field.  Once  the  flow  field  is  obtained  in  the  moving  frame  of 
reference,  the  results  are  compared  to  the  non-rotating  results 
to  exhibit  the  change  in  the  flow  field  appearance  from  the 
viewpoint  of  a  stationary  observer.  Of  particular  interest  are 
(1)  the  global  flow  structure  affected  by  the  rotating  motion, 
and  (2)  the  flow  field  properties  in  the  jet  plume. 


■:< 


SECTION  II 


GRID  REFINEMENT  STUDY 


1 .  ANALYSIS 

A  side-by-side  numerical  and  experimental  study  (Reference 
13)  of  the  entire  flow  field  for  a  configuration  consisting  of 
a  transverse  jet  issuing  from  a  blunt-nosed  ogive  cylinder  was 
recently  accomplished.  This  solution  will  be  designated  as 
case  A.  Figure  1  compares  the  schlieren  photograph  of  the  flow 
interaction  with  computed  density  contours.  The  global  fea¬ 
tures  of  the  flow  are  captured  by  the  computation,  although 
there  is  a  slight  discrepancy  in  the  projected  angle  of  the  bow 
shock  envelope  induced  by  the  jet.  Based  on  our  experience, 
this  is  a  first  clear  indication  of  insufficient  numerical 
resolution  for  this  complex  physical  phenomenon. 

Figure  2  depicts  the  comparison  of  the  experimental  and 
computed  pitot  pressure  profiles  at  three  locations  downstream 
of  the  jet.  The  shock  structure  can  be  determined,  but  as  is 
common  (References  13,  17,  23),  the  profile  is  smeared  due  to 
the  lack  of  grid-point  density  in  the  high-gradient  region.  We 
feel  t’  at  the  flow  field  solution  is  still  valid,  but  the  shock 
smearing  detracts  from  the  value  of  the  numerical  simulation 
for  this  comparative  study.  Another  poor  agreement  between 
computed  and  experimental  results  is  seen  near  the  wall, 
particularly  at  AX/D  =2.  We  are  reasonably  certain  that  some 
of  the  disagreement  is  due  to  insufficient  grid  resolution,  but 
this  observed  discrepancy  may  also  be  attributed  to  the 


combined  effects  of  the  inadequate  turbulence  modeling,  and 
possible  inaccuracies  in  the  experimental  data.  The  source  of 
error  in  pitot  pressure  measurements  might  be  traced  to  the 
fact  that  the  pitot  probe  was  aligned  only  in  the  freestream 
direction.  In  the  region  near  the  wall,  the  main  flow 
direction  is  dictated  by  the  jet  interaction  and  mixing.  The 
flow  deflection  to  the  unperturbed  flow  has  reached  a  value  as 

O 

high  as  45  .  This  particular  issue  of  the  data  accuracy  will 
remain  open  and  unresolved  in  the  present  analysis. 

As  a  final  representation  of  the  accuracy  of  the  computed 
solution,  we  examine  Figure  3,  comparing  the  computed  surface 
velocity  vectors  with  photographic  evidence  from  the  experi¬ 
ment.  Although  the  general  agreement  is  reasonable,  we  can  see 
a  major  topological  discrepancy  in  the  region  upstream  of  the 
jet.  The  experimental  oil  pattern  clearly  shows  two  lines  of 
convergence  and  two  lines  of  divergence.  The  computed  flow 
field  has  only  one  clear  line  of  convergence,  and  one  line  of 
divergence.  This  qualitative  disparity  represents  a  more 
serious  problem  in  numerically  replicating  three-dimensional, 
separated  flows  than  does  the  previously  described  quantitative 
di f  f  erence . 

On  the  basis  of  the  above  observations,  efforts  were  con¬ 
centrated  to  determine  the  requisite  grid  refinement  which  was 
sufficient  to  negate  the  discrepancies  noted  between  experi¬ 
mental  and  computational  data.  Particular  emphasis  was  placed 
on  the  resolution  of  the  separated  flow  upstream  of  the  jet  to 


achieve  a  grid  invariant  numerical  solution.  In  the  grid 
refinement  process,  we  constructed  the  grid  system  by  cluster¬ 
ing  grid  points  around  the  anticipated  resultant  shock  wave 
envelope  but  without  the  explicit  use  of  automated  grid 
adaptation . 

Four  grids,  designated  grids  A-D,  were  used  for  the  resolu¬ 
tion  study.  Side  views  of  each  grid  are  presented  in  Figure  4. 
All  grids  had  two  regions  of  clustering  in  the  normal  direc¬ 
tion.  Grid  points  were  clustered  near  the  body  to  resolve  the 
boundary  layer  and  away  from  the  body  to  capture  the  bow  shock 
wave.  The  most  coarse  grid  system,  used  in  the  previous  study 
to  calculate  case  A,  had  100  x  45  x  42  points  in  the  stream- 
wise,  normal,  and  azimuthal  directions,  respectively.  Grid  B, 
which  consisted  of  120  x  50  x  50  grid  points,  was  designed  to 
resolve  the  existence  and  location  of  the  possible  secondary 
vortical  structure  upstream  of  the  jet.  The  results  of  compu¬ 
tations  on  this  grid  are  designated  as  case  B.  The  additional 
20  streamwise  planes  are  added  in  the  region  upstream  of  the 
jet.  The  streamwise  spacing  in  this  region  is  decreased  by  a 
factor  of  two  over  that  of  grid  system  A  to  Order  6/10,  where  6 
is  the  undisturbed  boundary  layer  thickness  computed  at  the  jet 
location  on  the  opposite  side  of  the  cylinder.  The  third  grid, 
grid  C,  was  generated  in  an  attempt  to  capture  the  shock  loca¬ 
tion  downstream  of  the  jet;  it  had  110  x  70  x  50  points,  with 
the  majority  of  the  added  points  located  downstream  of  the  jet, 
between  the  body  and  the  experimentally  predicted  shock 
location . 


A  fourth  grid,  grid  D,  was  a  final  attempt  to  resolve  as 
many  of  the  features  of  the  flow  field  as  possible  in  one  grid 
system.  Grid  D  represented  increased  spatial  refinement  over 
cases  A,  B,  and  C,  both  upstream  and  downstream  of  the  jet. 

It  was  dimensioned  130  x  90  x  62,  and  was  fitted  to  the  general 
structure  of  the  predicted  shock  wave  envelope  of  jet  inter¬ 
action.  The  side  view  of  this  grid  system  shows  that  down¬ 
stream  of  the  jet,  the  grid  clustering  for  the  anticipated 
shock  wave  envelope  moves  away  from  the  body.  The  grid  expan¬ 
sion  is  limited  to  the  jet-injection  side  of  the  body;  the 
opposite  side  retains  the  same  basic  spacing  as  the  grids  in 
cases  A-C. 

Table  1  compares  cases  A-D  for  their  relative  grid  size, 
using  case  A  as  the  baseline.  Since  computational  time 
required  to  reach  convergence  is  proportional  to  the  number  of 
grid  points,  one  can  see  that  the  increase  in  resource  con¬ 
sumption  is  significant  for  all  three  grid  enhancements. 

The  numerical  solutions  were  obtained  using  the  three- 
dimensional,  mass-averaged  Navi er-Stokes  equations  with  a 

Table  1.  RELATIVE  GRID  SIZES 


Dimensions 


Relative  Size 


100 

x 

45 

X 

42 

120 

X 

50 

X 

50 

110 

X 

70 

X 

50 

130 

X 

90 

X 

62 

Relative  Size  = 


Number  of  Grid  Points/Number  of  Grid  Points 
for  Case  A. 


simple  flux-gradient  turbulence  model  (Reference  24).  The 
numerical  procedure  was  based  on  MacCormack ’ s  explicit 
predictor-corrector  algorithm  (Reference  25),  vectorized  for 

high  speed  processors.  The  data  processing  rate  was  4  x  10  ^ 
seconds  per  grid  point  per  iteration  on  the  Cray  2  at  NASA  Ames 
Research  Center.  Local  time  stepping  was  used  to  speed 
convergence  with  the  required  convergence  criterion  that  the 
consecutive  surface  pressure  change  less  than  one  percent  over 
a  characteristic  time  scale. 

2.  DISCUSSION  OF  RESULTS 

The  comparison  of  computed  pitot  pressure  profiles  with 
increasing  numerical  resolution  (cases  A,  C,  and  D)  together 
with  experimental  data  is  given  in  Figure  5.  In  general,  a 
steady  improvement  of  shock  wave  definition  was  gained  with 
higher  grid  point  density  and  better  distribution  of  mesh  nodes 
near  the  experimentally  predicted  location  of  the  shock.  In 
particular,  grid  system  D  represented  the  most  elaborate  grid 
refinement  and  yielded  the  best  agreement  with  data.  The  maxi¬ 
mum  deviation  between  experimental  data  and  numerical  results 
is  persistently  located  at  the  wave  front.  Further  refinement 
of  the  numerical  results  could  be  obtained  by  adding  more  grid 
points  and  by  redistributing  the  mesh  system  at  the  expense  of 
substantially  more  computer  resources.  The  use  of  shock  fit¬ 
ting  and/or  three-dimensional  adaptive  techniques  might  con¬ 
ceivably  be  even  more  effective  in  obtaining  the  desirable 
definition  of  the  shock  wave  envelope.  Finally,  the  influence 


of  the  turbulence  model  on  the  numerical  results  was  evaluated 
by  suppressing  the  numerical  transition  from  laminar  to  turbu¬ 
lent  flow  downstream  of  the  jet,  thereby  restricting  the 
numerical  simulation  to  laminar  flow  throughout  the  computa¬ 
tional  domain.  This  numerical  result  revealed  a  correct  out¬ 
ward  displacement  of  the  shock  wave  structure  downstream  of  the 
jet,  but  did  not  demonstrate  any  significant  change  in  the 
pitot  pressure  profile  near  the  body  surface. 

The  computed  surface  oil  flow  patterns  are  shown  in  Figure 
6  for  cases  A  and  B.  ine  two  numerical  results  are  markedly 
different  in  that  case  B  shows  a  distinct  secondary  convergence 
line  near  the  plane  of  symmetry,  whereas  no  such  line  is 
present  in  case  A.  The  locations  of  the  lines  of  convergence 
and  divergence  are  within  the  experimental  error-band  of 
surface  oil  streaks  under  the  identical  flow  conditions. 

A  side  view  of  the  vector  field  in  the  plane  of  symmetry 
ahead  of  the  jet  for  case  B  is  shown  in  Figure  7.  We  can 
clearly  detect  the  small  secondary  reversed  flow  region  adja¬ 
cent  to  the  body  surface.  The  greatest  dimension  associated 
with  the  secondary  structure  is  less  than  the  boundary  layer 
thickness  of  the  unperturbed  stream.  The  presence  of  the 
secondary  separation  region  affected  the  global  flow  field 
construct  negligibly.  This  can  be  seen  from  Figure  8,  which 
compares  the  plane  of  symmetry  density  contours  for  cases  A  and 
B.  One  may  conclude  that  the  increase  in  resolution  for  case  B 
does  increase  the  accuracy  of  the  description  of  the  separated 


flow  field.  The  predicted  secondary  vortical  structure  is 
found  to  replicate  experimental  observations.  The  secondary 
vortical  structure  influence  is  local,  so  that  the  relatively 
sparse  resolution  upstream  of  the  jet  is  not  overly  critical  in 
capturing  the  global  nature  of  the  flow  field. 

Figure  0  is  a  side  view  plot  of  the  velocity  vectors  in  the 
plane  of  symmetry  of  the  jet  for  case  D.  Compared  to  Figure  6, 
it  clearly  shows  that  the  solution  has  not  changed  in  this 
region  even  though  the  local  grid  refinement  has  been 
increased.  The  velocity  vectors  indicate  that  the  primary 
structure  contains  two  counter-rotating  vortices  upstream  of 
the  jet,  partitioned  by  the  primary  nodal  point  of  attachment. 
Downstream  of  the  nodal  point  and  immediately  upstream  of  the 
jet,  the  counter-clockwise  vortex  is  reinforced  by  the  jet 
stream.  Ahead  of  the  attachment  point,  the  fluid  which  origi¬ 
nated  far  upstream  and  immediately  adjacent  to  the  surface 
spirals  into  the  primary  clockwise  vortex  at  the  saddle  point. 
The  next  layer  of  fluid  similarly  forms  the  forward-most 
secondary  clockwise  system.  The  rest  of  the  entrained  fluid 
develops  into  the  secondary  counter-clockwise  rotating  vortex 
beneath  the  two  clockwise  vortices  and  appears  as  the  secondary 
separated  flow  region.  Figure  10,  which  presents  the  compari¬ 
son  for  the  overall  surface  shear  pattern  for  cases  B  and  D, 
clearly  shows  that  the  location  and  extent  of  the  secondary 
structure  are  identical  for  the  two  cases.  The  locations  of 
the  primary  convergence  lines  for  both  cases  agree  within  6/10, 


CASE  D 


Figure  9.  Velocity  Vectors  in  the  Plane  ot  Symmetry 


and  the  remaining  convergent  and  divergent  limiting  surface 
streamlines  show  even  closer  agreement.  The  overall  structure 
of  the  shear  pattern  remains  invariant  upstream  of  the  jet, 
with  no  new  topological  singularities  observed.  Thus,  for  this 
region,  the  solution  is  now  grid  independent,  and  the  solution 
of  case  B  stands  as  representative  of  the  complex  interaction 
structure . 


SECTION  III 


ROTATION  STUDY 

1 .  ANALYSIS 

When  the  equation  of  fluid  motion  is  expressed  in  the 
rotating  frame  of  reference,  the  fictitious  Coriolis  and 
centrifugal  forces  must  be  explicitly  included  in  the  formula¬ 
tion.  The  Rossby  number,  U/Lfl ,  a  convenient  measure  of  the 
ratio  of  inertial  and  Coriolis  forces,  becomes  an  additional 
dynamic  similarity  parameter  for  rotating  flow  systems 
(Reference  26).  For  a  typical  hypersonic  aerodynamic  config¬ 
uration  requiring  spinning  motion  to  maintain  stability,  the 
Rossby  number  is  around  750.  This  Rossby  number  was  simulated 
for  the  investigated  configuration  by  rotation  at  an  angular 
velocity  0  of  30  radians  per  second.  The  solution  obtained  at 
this  Rossby  number  will  be  termed  case  E.  In  order  to  deter¬ 
mine  the  influence  of  the  Coriolis  force  on  the  pattern  of  the 
flow  field,  solutions  were  also  obtained  at  Rossby  numbers  of 
75  (case  F)  and  7.5  (case  G) ,  spanning  a  range  of  three 
decades . 

Grid  A  was  chosen  for  the  rotation  study,  since  it  offered 
a  reasonably  accurate  solution  of  global  flow  field  features  at 
a  minimum  expenditure  of  computer  resources.  The  numerical 
algorithm  used  for  the  rotation  study  was  identical  to  the  one 
used  in  the  grid  refinement  study.  The  governing  equations  and 


the  numerical  boundary  conditions  were  altered  (as  follows)  to 
take  into  account  the  apparent  Coriolis  and  centrifugal  accel¬ 
erations  in  the  rotating  frame  of  reference.  The  Navier- 
Stokes  equations  can  be  written  in  the  flux  vector  form  as 
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When  the  above  equations  are  written  in  terms  of  the  rela¬ 


tive  velocity  in  the  rotating  frame  of  reference,  the  basic 
form  of  the  equations  remains  invariant.  The  additional 
acceleration  terms  in  the  rotating  frame  are  accounted  for  by 


adding  source  terms  to  four  of  the  five  equations.  This  is 
written  in  vector  form  as : 
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When  the  rotational  motion  is  restricted  about  the  streamwise 


axis  of  the  body  (ft  =  fti),  the  source  vector  has  only  three 
non-zero  terms : 

0 
0 

S  =  pCl7  y  +  2/oftw 
pQ?  z  -  2pftv 
pft  (vy  +  wz) 

It  may  be  interesting  to  note  that  the  Coriolis  force  has 
no  contribution  to  the  energy  equation,  since  it  is  a  force 
directed  at  right  angles  to  both  the  axis  of  rotation  and  the 
local  velocity  vector.  The  implementation  of  boundary  condi¬ 
tions  at  the  solid  surface  must  also  be  changed  in  the  rotating 
frame.  The  compatibility  condition  at  the  surface  can  be 

satisfied  by  changing  the  normal  APwa^  =  0  to  APwaj^ 

where  p  and  r  are  evaluated  at  the  body  surface. 

The  relative  impact  of  the  Coriolis  and  centrifugal  forces 
can  be  determined  by  computing  and  comparing  their  magnitudes 
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throughout  the  flow  field.  The  centrifugal  forces  are  linearly 
proportional  to  the  radial  distance  from  the  axis  of  rotation, 
and  their  effect  should  become  stronger  in  the  outer  region  of 
the  flow.  The  Coriol  is  forces,  however,  must  reside  where 
there  are  high  velocity  components  normal  to  the  axis  of  rota¬ 
tion.  It  follows  that  the  strongest  Coriolis  force  should 
occur  as  the  jet  issues  from  the  body.  Any  small  change  in 
structure  at  this  point  should  theoretically  propagate  and 
amplify  downstream  in  the  trajectory  of  the  jet,  so  the  jet 
trajectory  may  serve  as  the  best  indicator  of  the  presence  of 
significant  Coriolis  effects. 

In  the  interaction  region  immediately  upstream  and  down¬ 
stream  of  the  jet,  velocities  are  small  compared  to  other 
regions  of  the  flow  field,  especially  near  the  surface.  The 
separated  flow  pattern  is  known  to  be  sensitive  to  any  small 
disturbances.  We  could  expect  that  rotational  forces  might 
cause  significant  changes  in  the  local  structure.  We  will 
examine  the  local  shear  stress  formations  upstream  and  down¬ 
stream  of  the  jet  for  possible  asymmetries. 

2.  DISCUSSION  OF  RESULTS 

Figure  11  depicts  a  comparison  of  the  relative  magnitude  of 
the  maximum  Coriolis  and  centrifugal  forces  for  cases  E,  F,  and 
G.  The  centrifugal  force  is  three  orders  of  magnitude  weaker 
than  the  Coriolis  force  for  case  E,  the  lowest  spinning  rate 
simulated.  Even  at  an  angular  velocity  100  times  greater  than 
the  case  of  interest,  the  centrifugal  force,  which  is  propor- 
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tional  to  fi  ,  is  an  order  of  magnitude  smaller  than  the 
Coriolis  force.  For  all  three  cases  studied,  the  centrifugal 
force  probably  has  not  substantially  affected  the  formation  of 
the  f low  field . 

Figure  12  consists  of  two  contour  plots  which  show  the 
magnitude  of  Coriolis  force  for  sections  of  the  entire  flow 
field  for  case  E.  In  the  longitudinal  plane  of  symmetry,  the 
Coriolis  forces  are  shown  to  be  strongest  at  the  jet  exit,  and 
negligible  elsewhere.  The  contour  levels  exhibited  in  the 
figure  represent  the  upper  90  percent  values  of  the  magnitudes 
of  Coriolis  force  for  the  complete  flow  field.  The  cross- 
sectional  contours  taken  at  the  jet  location  again  show  that 
the  Coriolis  forces  are  only  significant  near  the  jet  exit. 

Figure  13  compares  the  cross  sectional  view  of  the  Mach 
number  contours  at  the  jet  location  for  cases  A,  E,  F,  and  G. 

We  see  that  case  E,  which  satisfies  the  Rossby  number 
similarity  for  a  typical  hypersonic  configuration  of  practical 
interest,  is  identical  to  case  A.  Rotation  at  a  rate  10  times 
higher  than  case  E  still  yields  no  appreciable  difference.  It 
is  only  when  the  angular  velocity  is  increased  by  two  orders  of 
magnitude  over  the  rate  of  interest  that  a  noticeable  change  in 
the  shock  structure  takes  place.  In  this  case,  the  shock  enve¬ 
lope  immediately  acquires  an  asymmetrical  shape.  The  trajec¬ 
tory  of  the  jet  is  disturbed  only  in  case  G,  where  it  already 
has  begun  to  shift  away  from  the  direction  of  rotation  due  to 
the  restoring  effect  of  the  Coriolis  force. 
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Downstream  of  the  jet,  the  effects  due  to  rotation  are 
consistent  with  previous  findings.  Since  the  most  pronounced 
rotational  effects  occur  at  the  jet  exit,  the  jet  plume  should 
accentuate  disturbance.  Figures  14  and  15,  which  present  Mach 
contours  at  stations  six  and  twelve  jet  diameters  downstream  of 
the  jet,  support  this  observation.  At  the  rotation  rate  that 
satisfies  the  similarity  condition,  no  rotational  effect  can  be 
observed  on  the  flow  field  Only  at  the  highest  rotation  rate 
computed  is  there  an  appreciable  change  in  the  shock  structure 
and  jet  trajectory.  For  this  case,  the  disturbance  originates 
at  the  jet  exit  and  propagates  downstream. 

The  behavior  of  the  shear  vectors  is  identical  for  all  four 
cases  in  the  near  surface  region  upstream  of  the  jet.  This  is 
not  surprising  in  view  of  the  computed  magnitude  of  CorioliB 
forces  outlined  earlier  In  the  separated  flow  region  upstream 
of  the  jet,  there  are  no  rotat  1  on  forces  with  sufficient 
strength  to  alter  the  flow 

Results  are  slightly  different  downstream  of  the  jet, 
however,  as  shown  in  Figure  16  The  surface  shear  patterns  are 
the  same  for  cases  A,  E,  and  F' ,  but  case  G  begins  to  show  some 
asymmetry  due  to  rotation  Tie  footprint  of  the  jet  wake  is 
shifted  slightly  away  from  the  direction  of  rotation,  consis¬ 
tent  with  the  shift  in  both  the  jet  trajectory  and  the  shock 
structure  . 

The  existence  of  significant  changes  in  the  flow  properties 
near  the  jet  can  be  determined  by  examining  the  close  up 
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behavior  of  fluid  density  in  that  region.  Figure  17  shows  the 
cross-sectional  contour  of  density  at  the  jet  location  in 
detail  for  cases  A,  E,  F,  and  G.  The  contours  for  cases  A,  E, 
and  F  are  identical  within  the  range  of  numerical  error.  Case 
G,  simulating  the  highest  rate  of  rotation,  shows  large  changes 
of  the  local  flow  field  properties  induced  by  the  much  higher 
rotation  rate.  The  asymmetry  of  density  contours  in  case  G 
indicates  a  changes  in  trajectory  of  the  jet  due  to  the 
Coriolis  forces.  Note,  however,  that  very  near  the  jet  exit 
the  density  contours  retain  the  symmetric  distribution. 

Figure  18  shows  the  side  view  of  the  plane  of  symmetry  of 
the  jet  for  cases  A,  E,  F,  and  G.  Here,  no  changes  are  noted 
in  the  flow  properties  near  the  jet  for  any  of  the  computed 
results.  This  observation  was  anticipated  and  used  as  a 
benchmark  for  the  present  investigation.  Since  the  centrifugal 
force  is  negligible  in  comparison  with  the  Coriolis  force  and 
the  latter  must  be  directed  at  right  angles  to  both  the  axis  of 
rotation  and  the  jet  velocity,  the  rotational  forces  must  be 
exerted  only  in  the  circumferential  direction.  The  effect  of 
rotation  projected  in  the  plane  of  symmetry  of  the  jet  is 
neither  expected  nor  generated  numerically  for  the  cases  under 
consideration . 
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SECTION  IV 


CONCLUSIONS 

Flow  fields  of  a  supersonic  jet  ejected  from  an  ogive 
cylinder  into  a  hypersonic  freestream  with  and  without  rotating 
motion  have  been  numerically  simulated.  For  the  non-rotating 
calculations,  the  three-dimensional  separated  flow,  consisting 
of  two  pairs  of  counter-rotating  vortices,  replicates  experi¬ 
mental  observations.  Even  though  the  definition  of  the  resul¬ 
tant  shock  wave  is  smeared  by  the  shock  capturing  scheme,  the 
computed  complex  vortical  structure  upstream  of  the  jet  is 
accurate  and  independent  of  further  mesh  system  refinement. 

The  rotating  flow  system  was  computed  in  a  rotating  frame 
of  reference.  Three  different  body  rotation  rates  correspond¬ 
ing  to  a  range  of  Rossby  numbers  of  750,  75,  and  7.5  were 
examined.  At  the  highest  Rossby  number  corresponding  to  the 
dynamic  condition  of  practical  engineering  interest,  the  flow 
field  in  the  rotating  frame  of  reference  is  identical  to  that 
in  the  non-rotating  case.  The  entire  flow  formation  rotates  as 
a  solid  bulk.  The  jet  plume  and  the  downstream  shock  wave 
envelope  exhibit  a  circumferential  shift  due  to  the  body 
rotation  only  at  the  lowest  Rossby  number  examined,  7.5. 
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LIST  OF  SYMBOLS 


D  diameter  of  jet  aperture 

e  total  energy  per  unit  mass 

F,G,H  flux  vectors  defined  by  equation  1 
k  thermal  conductivity 

L  characteristic  length 

M  Mach  number 


pressure 


pitot  tube  pressure 


r  radial  distance  from  the  axis  of  symmetry 

R  radius  of  cylinder 

Rey  Reynolds  number 

5  rotational  source  vector  defined  by  equation  2 

t  time 

T  temperature 

u,v,w  velocity  components  in  the  Cartesian  frame 
velocity  vector 

x,y,z  coordinates  in  the  Cartesian  frame 

6  boundary  layer  thickness 

p  fluid  density 


stress  tensor 


angular  velocity,  radians  per  second 


Subscript 


i,j,k  unit  vector  in  the  Cartesian  frame 


